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Given a standard Brownian motion _B'' = (B^)o<t<T with drift /i ^ IR and letting g denote 
tile last zero of Bf^ before T, we consider the optimal prediction problem 

14= inf Elg-rl 
0<T<T ' ' 

where the infimum is taken over all stopping times r of . Reducing the optimal prediction 
problem to a parabolic free-boundary problem and making use of local time-space calculus 
techniques, we show that the following stopping time is optimal; 

T. = inf {t G [0,T] I Bj'' < b-{t) or B^' > b+{t) } 

where the function t i— > 6_ (t) is continuous and increasing on [0, T] with f)_ (T) = 0, the 
function t ^ b^{t) is continuous and decreasing on [0,T] with f)-i-(T) = 0, and the pair fe_ 
and fe+ can be characterised as the unique solution to a coupled system of nonlinear Volterra 
integral equations. This also yields an explicit formula for in terms of b_ and b_|_. If = 
then fe_ = — fe+ and there is a closed form expression for b± as shown in lid using the method 
of time change from 0|. The latter method cannot be extended to the case when ^ and 
the present paper settles the remaining cases using a different approach. 
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1. Introduction 

Imagine a stochastic process X = {Xt)o<t<T observed over the time interval 
[0, T], and let g denote the last zero of X before T. Formally, we set g = 
sup{t G [0,T] I Xt = 0}so that Xg = when X is left-continuous or Xg^ = 
when X has left limits (and the set is non-empty in either of the two cases). 
Clearly, at any time t < T the value of g is unknown (unless X is trivial), and 
it is only at the terminal time T itself that we know when the last zero g oi X 
occurred. However, this is often too late: typically one wants to know how close X 
is to g at any time t < T and then take some action based on this information. 
In the present paper we take up this question when X is a standard Brownian 
motion = (i?Oo<t<i with drift fi £ M and search for a stopping time r^, of B^^ 
that is as 'close' as possible to g. The optimal prediction problem is more precisely 
formulated in ()2.2p below and the main result is given in Theorem 1. 
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It was recently observed in [11] that the optimal stopping time in (2.2) when /u = 
is equal to = { t S [0, T] | 15^1 > z^\jT—i } where > is a specified constant 
and B = {Bt)o<t<T is a standard Brownian motion (with no drift). This stopping 
time is equal in law to the stopping time f^: = inf {t^[0,T]\St — Bt> z^:\jT—t } 
which is known (cf. [4]) to minimise ^{St — -Br)^ over all stopping times t oi B 
taking values in [0, T], where we set St = maxo<s<t Bg for t € [0, T]. Therefore when 
there is no drift, the space domain problem of stopping as close as possible to the 
ultimate maximum of Brownian motion, and the time domain problem of stopping 
as close as possible to the last zero of Brownian motion, are exactly equivalent. It 
is natural to consider whether this connection is preserved when the drift fi ^ 0. 
In this case it is known (cf. [3]) that the solution to the space domain problem (of 
minimising ^{St^-Bi^f) is given by = inf {t G [0,r] | hi{t) < S^-B^ < b2{t)} 
where bi and 62 are specified functions of time and G [0,T). This stopping time 
has a distribution which is rather different from that of above (when /j, is away 
from zero) raising the question of whether the solution to (2.2) below will bare any 
resemblance to it. In this paper we show that the optimal stopping time in (2.2) is 
in fact more similar to r* above. Predicting the last zero of Brownian motion with 
drift is therefore quite distinct from predicting its ultimate maximum. 

Turning to a general stochastic process X = {Xt)o<t<T one can ask similar 
questions about predicting the last zero of X in the interval [0,T]. Certainly, if 
X is discontinuous or multidimensional one would expect the problem to be much 
harder. Although these issues are not directly addressed in the paper, a simple 
approach based on the local time of X is briefly presented in Remark 1 below. 
This may be used as a starting point from which to study similar problems for 
more general processes. 



2. Formulation of the problem 

We begin our exposition by formally introducing the setting and the problem to be 
studied. Let B = {Bt)o<t<T be a standard Brownian motion defined on the proba- 
bility space (0, J^, P) with Bq = under P. For any /i G iR let B^' = (Sf)o<f<r be 
the Brownian motion with drift fi given by B^ = jit + Bt for t £ [0, 1] . Define the 
last zero of B^ as 

g = sup{0<t<T\B^ = 0} (2.1) 

so that g is well-defined under P since S'^ starts at zero under this measure. Later 
on we shall need to allow B^ to start at arbitrary values in M, raising the question 
of which value to assign to g if B^ does not hit zero at all on [0, T]. However, as it 
turns out (see Lemma 1 below), the value of g in this case is irrelevant. 
Consider the optimal prediction problem 

K= inf E|g-r| (2.2) 

o<T<r 

where the infimum is taken over all stopping times r of B^^. The gain process 
{\9 ~t\)o<t<T is not adapted to the natural filtration {J-'f^)o<t<T of B^^ as g is 
only J^^'^ measurable, meaning that the problem above falls outside the scope of 
standard optimal stopping theory. However, it is well known that there exist close 
links between g and the maximum of S^, and these we now exploit with the aim 
of reducing (|2.2p to an equivalent optimisation problem to which the standard 
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techniques of optimal stopping for Markov processes (see e.g. [7|]) can be applied. 

Let = {Sf)o<t<T be defined by = maxo<s<t-B^ so that S^" denotes the 
maximum of B^' up to time t G [0,r]. Then it follows (cf. & [1]) that the 
distribution of is given by the formula 



X — fj,t 



-X — fit 

"7^ 



(2.3) 



for all {t,x) € [0, oo) x where ^{x) = J^^ip{z) dz denotes the distribution 

function of a standard normal random variable and (p{x) = (l/\/27r) e~^^/^ denotes 
its density for x £ M. 

Lemma 2.1: The optimal prediction problem (|2.2|) above is equivalent to the 
standard optimal stopping problem 



V= inf E 

0<r<T 



H{t, B^) dt 



(2.4) 



where the function H is given by 

H{t,x) = 2F'^~^\T-t,x)I{x > 0) + 2F^^\T-t,-x)I{x < 0) - 1 (2.5) 



for all {t,x) G [0,T]x]R, and where 14 from 



is given by = V + E{g)- 



Proof : We follow the same approach as in Lemma 1] and [9|] . Fix any stopping 
time T of taking values in [0,T]. We then have 



\g -t\ = {t - gy + {t - g) = {t - g)+ + g -T Ag 

f'T rr 

= I{g<t) dt + g- I{g> t) dt 
Jo Jo 

= g+ I [21{g<t)-l]dt. 



(2.6) 



From Fubini's theorem one finds that 



/ I{g<t)dt = I E[I{T>t)I{g<t)]dt 
Jo . Jo 



(2.7) 



/ E[/(r>t) 
Jo 



I{g<t)\Tndt 



P(5 < tl^Drft 



Using the continuity of the sample paths of B^ to compute the conditional proba- 
bility above we arrive at the following expression: 



P( min > I J^f ") if Bl > 



P{g<t\ Tt) 



t<r<T 



P( max B^^ <0\ J^r) if < 



(21 



for t £ [0, T] (the case B[^ = can be ignored since the amount of time Brownian 
motion spends at zero has Lebesgue measure zero). If we consider the case when 
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-Bf > we find that 



P{mmB^ >0\:Fr) = P( min (i^A^ _ ^jM) + > 0\ J^^) (2- 



P( max (-B^ + Bf) <b\n 



b=Bf 



for any t £ [0, T] since i?^ has stationary independent increments. A similar argu- 
ment when B^ < shows that 



P( max i^A. < I = P(5^_, < -b) 



b=B'f 



(2.10) 



for any t € [0, T]. Taking expectations on both sides of (|2.6p and using (|2.7p - ()2.10p 
we get 



E|5-r| = E(<7) + E 



2F^'^\T-s,B^)I{B^ > 0) 

+ 2F(^)(r-s, -B^) I{B^ < 0) - 1 



(2.11) 



ds 



for any stopping time r of B^^ with values in [0,r]. Taking the infimum over all 
such stopping times we conclude the proof. □ 

In order to apply the standard techniques from the theory of optimal stopping 
for Markov processes (see e.g. 0]) it is necessary to extend the problem ()2.4p by 
allowing B^^ to start at any time t E [0, T] at any point x in the state space. 
Since B^^ is a time-homogeneous Markov process there exists a family of Markov 
measures {Px)x€iR such that Px(i?o = x) = 1 for any x G M. Using these measures 
to change the starting point of the process results in the problem (|2.4p becoming 



V{t,x) 



inf I 

0<T<T-t 



H{t+s,B^)ds 



(2.12) 



for any (t, x) £ [0, T]x M so that the value 14 of the optimal stopping problem 
()2.2p is given by = ^(0,0) + E(g). However, it is well known that the process 
{x+B^)t>o under Pq provides an explicit realisation of B^^ under P^; for any x G M. 



This together with the fact that the optimal stopping time in (j2.4p is the first 
hitting time to a set (this follows from general theory and will also be proved 
below) leads to the further simplification 



V(t,x)= inf E 
o<T<r-t 



[ H{t+s,x+B^)ds 
Jo 



(2.13) 



where P is the original measure Pq under which B^ starts at zero at time 0. This 
expression then is the form we shall use when solving the optimal prediction prob- 
lem (|2.2p . Note that V < since we can always insert the stopping time identically 
equal to zero in the expectation above. 
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3. Result and proof 

In order to state our main result we need the following definitions. Let the function 
/ denote the density of B^" so that 

/(t,x) = -i=e-^(--'^*)^ (3.1) 

V ZTTl 

for all {t,x) G (0,T] xM. Recalling H in (|2.5p above define the function K as 

K{t,x,s,z^,z+) = E[H{t+s,x+B^)I{z^ < x+B^ < z+)] (3.2) 

^ ^ H{t+s,x+b)f{s,b) db 

X+Z- 

for all {t, x) £ [0, T] x M, all s £ [0,T -t] and z_ < z+ in M. Lastly the 
set {H < 0} := { {t,x) G [0,r] xiR | H{t,x) < 0} wih play a prominent role in 
our discussion. A direct exam- ination of H reveals the existence of two con- 
tinuous functions /i_ and mapping [0, T] into M such that {H < 0} = 
{{t,x) S [0,T]x M \ h^[t) < X < h^{t)} . Furthermore is increasing and /i+ is 
decreasing with /i_ < < /i+ and h^[T) = h^{T) = 0. We may now state our 
main result. 

Theorem 3.1 : Consider the optimal stopping problem (I2.13p . Then there exist 
continuous functions 6_ and b^ on [0, T] such that the optimal stopping set is given 
by 

D = {{t,x) £ [0,T]xlR \ X < b^{t) or x> b+{t) } (3.3) 

where t b^{t) is increasing, t ^ b^{t) is decreasing, and b-(T) = 6+(T) = 0. 
This means that the stopping time 

TD{t, x) = inf { s e [0, T-t] I x+Bi^ £ D } (3.4) 

is optimal for all {t,x) £ [0, T] xM. The value function V from ()2.13p is given by 

l-T-t 

V{t,x)= K{t,x,s,b-{t+s),b+{t+s)) ds (3.5) 
Jo 

for all (t, x) £ [0, T] x M, and the functions 6_ and b^ themselves are uniquely 
characterised by the coupled system of non-linear Volterra integral equations 

.T-t 

/ K{t,b^{t),s,b^{t+s),b+{t+s))ds = Q (3.6) 
Jo 



rT-t 

/ K{t,b+{t),s,b^{t+s),b+{t+s)) ds = 
Jo 



(3.7) 



in the sense that they are the unique solution to (j3.6p + (j3.7p in the class of con- 
tinuous functions t b-{t) and t i— > 6+(i) on [0,T] satisfying b-(t) < h^(t) 
and &+(t) > hj^{t) for all t £ [0, T]. Finally, the value from (j2.2p is given by 
V-t = y(0, 0) -|- E(gf) and the optimal stopping time for this problem is r/)(0, 0). 
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Figure 1. A computer drawing of the optimal stopping boundaries for Brownian motion with drift 
m < and fi2 > 0. The dotted line is the optimal stopping boundary for Brownian motion with zero 

drift. 



Proof : 1. Existence of optimal stopping time. We begin by showing that an op- 
timal stopping time for the problem (j2.13p exists. Since H is continuous and 
bounded, and the flow x x + is continuous, it follows that the map 
(t, x) I— > H{t + r,x + Br) dr] is continuous and thus upper semicontinuous 

(use) for any stopping time r with values in [0,T]. The function V is therefore use 
as well (being the infimum of use functions) and so by general results of optimal 
stopping (see Corollary 2.9]) it follows that an optimal stopping time for the 
problem (j2.13p exists. Moreover, this stopping time is given by (j3.4p above where 
the stopping set equals L> = { (t, x) G [0, T\x]R \ V{t, x) = } and the continua- 
tion set is given by C = { (t, x) G [0, T] x iR | V{t, x) < } =0". The fact that D 
is closed (and C is open) follows from the fact that V is use. 



2. Shape of D. From the functional form of H given at ()2.5p (see also (12. Sp 
above), one sees that the map 1 1-^ H{t,x) is increasing on [0,r] for every x ^ M. 
Fix any x £ IR and s < t in [0, T] and set Tg = td{s,x) and = T£){t,x). Since 
< Tf < T—t < T — s one finds that 



V{t,x) - V{s,x) 



(3.8) 



> E 



H{t+r,x + B^)dr 



H{s+r,x + B^)dr 



H{t+r, x+B^!-) - H{s+r,x+Bi;) dr 



> 



from where we derive the important fact that 

1 1— > V{t, x) is increasing on [0, T] 



(3.9) 



for each x £ M given and fixed. A direct consequence of this is that if any point 
(t, x) e D, then all the points {t+s,x) £ D for s £ [0,T—t] since > V{t+s,x) > 
Vit,x) = 0. 
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Turning to the continuation set C, take any point (t, x) in the open set {H < 0} 
and let U C {H <0} be any open neighbourhood of {t,x). Defining the stopping 
time ajj as the first exit time from U, we see that 



V{t,x) < E 



H{t+s,x+B^)ds 



< 



(3.10) 



showing that {H <0} C C. We now show that D has the form given in (|3.3p . Fix 
any t in [0, T) and x > and suppose that (t, x) G D. For any y > x consider the 
process B^^ started at {t,y). Since {s,x) £ D for all t < s < T it fohows that Bf^ 
started at (t, y) must enter D before crossing the level x, in other words y + Bg > x 
for all < s < TD(t,y). From equation ()2.5p we see that the map x H{t,x) is 
increasing on iR_|_ for each t £ [0, T] given and fixed. Hence for y > x > we have 



V{t,y) = E 



> E 



rn{t,y) 



H{t+s,y+B^)ds 



rD{t,y) 



H{t+s,x + B^)ds 



(3.11) 



> E 



Tn{t,x) 



H{t+s,x+B^)ds 



since Toit, x) = 0. Therefore if {t, x) £ D with x > 0, then (s, y) £ D for all y > x 
and s £ [t, T]. For x < a similar argument can be made showing that if (t, x) £ D 
then {s,y) £ D for all y < x and s £ [t,T]. It follows then that D has the form 
given in (13. 3p where the functions 6- and 6+ are defined as 



6+(t) = inf {x > I (t,x) e -D} and 
6_(t) = sup{x < I {t,x) £ D} 



(3.12) 



for all t £ [0, T]. We also see that 6+ is decreasing and 6_ is increasing. 

3. Functions 6+ and 6_ are finite valued. We show that the optimal stopping 
boundaries 6_|_ and b_ are finite valued, which also shows that the stopping region 
D is strictly greater than the set { (T, x) \ x £ M} . We shall only present the proof 
for 6+ since the argument for 6_ is analogous. 

Suppose that the function b-^-{t) is not finite valued for all t £ [0, T] and define the 
time t^ £ [0, T] as = sup{t £ [0, T] | 6+(t) = oo}. There are two possibilities: 
either 6+ has a jump discontinuity at t^, jumping down from infinity to a finite 
value (including the case when t^^ = T), or 6+ has an asymptote at i,,. Suppose 
first that 6-1- jumps from infinity at t^ and that £ (0,T). If Tx = T£){0,x) and 



ax = inf {t £ [0, T] | x+B^ < /i_|_(0) + l } , then it is clear that 
as x ^ oo. Therefore we see that 



and 



T 



> V{0,x) 



H{s,x+B^)ds I{Tx<ax) 



+ E 



/' 

Jo 



H{s,x+B^)ds I{Tx>ax) 



(3.13) 



> E[eTxI{Tx<ax)] - TE[l{Tx>ax)] 

where e := ml {H{t,x) \ (t, x) £ [0, T] x [/i-|-(0) + l, cxd) } > and where we have 
used the fact that H > —1. Letting x ^ oo and using the dominated convergence 
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theorem we obtain 



> lim V{0,x) >eU>0 



(3.14) 



which is a contradiction. This shows that we cannot have t* E (0,7"). However, 
since is right-continuous (this will be shown below) we cannot have t^, = 
either, and so 6+ cannot have a jump discontinuity at £ [0, T). Finally, if = T 
then letting a; — > cxd and using the dominated convergence theorem we find that 



> V{0, x) = lim E 



rT 



H{s,x+B^)ds 



(3.15) 



[ H{s,oo)ds = T > 
Jo 



which is a contradiction. This excludes the case t^, = T and establishes the claim. 

Suppose now that asymptotes to infinity at and that G (0,7"). It then 
follows that < liminfi.^oo "^z < limsup^.^o^ r^. < + 5 < T for some 5 > 0, 
and the same argument as in (j3.13p yields > limsup^^o^ y(0, a;) > et* > 0. 
Lastly, suppose that 6+ asymptotes to infinity at 0. By extending the terminal time 
to r' > T and considering our problem on the interval [0, T'] instead of [0, T] we 
will (since 6+ is decreasing) shift to be strictly positive, reducing it to the case 
already considered. Therefore 6+(t) must be finite for all t € [0, T]. 

4. Continuity of V. We show that (t, x) i— > V{t, x) is continuous on [0, T] x ]R. 
For this, take any t G [0,r] and x,y £ M, set = T£){t,x) and Ty = Toit^y) and 
suppose without loss of generality that x < y. It follows then that 



H{t+s, y+B^) - H{t+s, x+B^) ds 



<V{t,y)-V{t,x) < E 



H{t+s, y+B^) - H{t+s,x+B^) ds 



(3.16) 



Turning to equation (|2.5p and manipulating the quantities we find that 



Hit,y)-Hit,x) 

= 2(^F(-^)(r-t,y) -F(-'^)(r-t,2:)) I{x > 0) 
- 2(F('^)(r-t, -x) - F^^^\T-t, -y)) I{y < 0) 
+ 2( F^-^'\T -t,y) - F^^'\T -t, -x)) l{x<0<y). 



(3.17) 



Applying the mean value theorem in the first and second line upon noting that 



4^\t,x) 
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for all {t,x) S (0,T] xiR+ we obtain the inequality 



1 



4:(y-x) . , 

<H{t,y)-H{t,x) 



+ 1^1 I /(y < 0) - 21{x<0< y) 



(3.19) 



+ /(x > 0) + 21{x <0<y) 



which we can insert into equation (|3.16|) to get 



\Vit,y)-Vit,x) 



(3.20) 



< E 



T-t 



4(y-x)((r-t-.)-i/2 + |;x|) 
+ 2 1{x+B^ < < y+B^) ) ds 
4(y-x) (2VT^t + \fi\{T-t) 



< 4(y-x) (2VT^t + |/i|(r-t)j + (y-x) 

< (uVf+iT\f,\){y-x) 



upon using that ^(z) — ^{w) < z — w for w < z in JR. Hence we conclude that 
x I— > V{t, x) is continuous on M uniformly over all t G [0, T]. 

It remains to show that t V{t, x) is continuous on [0, T] for every x G M. To 
see this fix x S iR, take any s < t in [0, T] and set = T£)(s, x) and Tt = T£)(t, x). 
Using the stopping time a = Tg A (T—t) instead of T£){t,x) one finds that 



< V{t,x) - V{s,x) 
<E [ H{t+r,x+B^!) - H{s+r,x+B^^)dr 



(3.21) 



H{s+r,x + B'')dr 



Taking the limit as t — s ^ and using the dominated convergence theorem, the 
continuity of t i— > H{t, x) and the fact that < Tg — a < t — s, we see that both 
expectations on the right hand side vanish and we conclude that (t, x) i— > V{t, x) 
is continuous on [0, T] x JR. 

5. Free-boundary problem. We now formulate a free-boundary problem that the 
value function V solves. This differential equation will be useful to us later on, 
but is also interesting in its own right and can be used as the departure point in 
computing numerical values for the boundaries 6+ and 6_ and for the value function 
V. It is well known from the theory of Markov processes (see e.g. 0, Chapter III, 
Section 7]) that V is C^'^ in the continuation region C and satisfies the following 
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version of the Kolmogorov backward equation: 

Vt{t, x) + x) + ^V^x{t, x) = -H{t, x) for all {t, x) e C (3.22) 

Vit, x) = for all (t, x) £ D. (3.23) 

The last condition which forms part of the free-boundary problem is the smooth 
fit property which reads 

X I— > Vx{t,x) is continuous over b-{t) and fo+(t) for all t € [0,T). (3.24) 

This result, which we establish below, will play an important role in the derivation 
of the integral equations (j3.5p and (j3.6p + (|3.7p . 

6. Smooth fit. We show that x i— >• V^(i, x) exists and is continuous over the optimal 
stopping boundaries ?)_ and Since the arguments for the two curves are very 
similar we will only present the proof for . 

Fix any time t G [0,T) and e > 0, and set x = 6+(t) and = T£){t, x — e). From 
equations (j3.16p and (|3.19p we find that 



< V{t,x) - V{t,x-e) 
< E 



(3.25) 







H{t,x+B^) - H{t,x-e+B^)ds 
4e((r-f-s)-i/2 ^ )+2I{0< x+B^ < e) 







ds 



2{VT^- ^T-t-Ts) +\fJ,\Te 

1(0 < x+B^ < e) ds 



< E 

< 4eE 
+ 2E 



Consider the final expectation above and denote by i^{B'^) = (£f (i?'^))^^^ the local 
time of B^ at the level x. Then it follows from the occupation times formula (see 
e.g. [1, p. 224]) and the integral mean value theorem that 



1 

e Jo 



/(O < < e) ds 



e Jo 



I{-x <a< -x+e) i'^^ {B^) da (3.26) 
:^+^{^B^')da = l-^+'^^{B^') 



where jy^ is a random variable taking values in [0,e]. Define the stopping time 
to be the first time i?^ rises above the level x when B^^ starts at time t at the 
point x—e, in other words Ue = inf { s G [0,T— t] | x—e+Bs = x} . Since the map 
t ^ b^{t) is decreasing it follows that < < (T^. Recalling that t i-^ —fit is a 
lower function of Brownian motion at 0+ we see that as , and consequently Ts , tend 
to zero as e 1 0. 

Dividing equation ()3.25p through by e and passing to the limit as e | 0, we 
see that the first expectation disappears (upon using the dominated convergence 
theorem) so that 



0<lim^(^'-)-^^^'---^)<2 1imE 
£^0 £ e-*0 



(3.27) 
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where for the final equahty we use the fact that l^^~^^'[B^) is bounded by an 
integrable random variable when e > and the well-known result that (t, a) i— > 
£f{B'^) is continuous on x M. The former can be seen by recalling the Tanaka 
formula 



i^r-ai 



\a\ + [ sign(5^-a)dS^ + ^^(S^) 
Jo 

\a\+Pt + f^ [ sign{B^ -a) ds + e'i{B^) 
Jo 



(3.28) 



where a G iR is any number and Pt = /q sign(i?^ — a) dBg for t > is a standard 
Brownian motion by Levy's characterisation theorem. It follows from (j3.28p that 



sup sup4"(S^) < max \Bt\ + max \f3t\ +2\fj.\T 

0<t<TaeIR 0<t<T 0<t<T 



(3.29) 



and the claim above follows by noting that the right-hand side in (|3.29p defines an 
integrable random variable. We therefore conclude that x V{t, x) is differentiable 
at b+{t) and that Vx{t,b+{t)) = 0. 

A small modification of the argument above shows that x Vx{t,x) is contin- 
uously differentiable at b^{t). Indeed, taking 6 > and setting ts = r£)(t, x — (5) 
where x = b+ (t) , we find using equations (|3.16p , (|3.19p and (|3.26p that 



V{t,x-6+e) - V{t,x-5) 



< 4E 



2{VT^ - y^T-t-n ) + Ts] (3.30) 



2E 



where e G (0, 6) is a given number and % is a random variable taking values in 
[0, e]. In exactly the same way we find that 



4E 



2{VT^t-y'T-t-Ts)+\ij\ 



2E 



(3.31) 



V{t,x-5) - V{t,x-5-e) 



where e > is a given number and pe is a random variable taking values in [0,e]. 
Letting first e J, in (|3.30p and (j3.3ip (upon recalling that V is C^'^ in C so that 
Vx{t,x — 5) exists) and then (5 | in the resulting inequalities (upon recalling that 
TS 0) we see that the map x i-^ Vx{t,x) is continuous at 6+(t). 

7. Continuity of boundaries. We show that the functions t b-{t) and 1 1— > 6+(t) 
are continuous on [0, T]. As before we will restrict ourselves to 6+ since the argument 
for b- is similar. 

We first show that b^ is right-continuous. Fix any t S [0,T), let tn I t and 
consider the limit 6+(t+) := lim„_>oo b+{tn) which exists as 6+ is decreasing. Since 
(tn,b+{tn)) G D for all n > 1 and D is closed, it follows that (t,6+(t+)) € D and 
so from (|3.12p we see that 6+(t+) > 6+(t). On the other hand, the fact that 6+ is 
decreasing implies that &+(t) > 6+(t„) for all n > 1, and passing to the limit as 
n ^ oo we obtain the reverse inequality. 

We now show that 6+ is left continuous. For this, suppose there exists some 
t G (0,T] at which b+{t—) > b+{t) and choose any x G (6+(t), ^^(t— )) . Since 
b+ > /i+ and is continuous, it follows that x > /i+(s) for all s £ [si,t] for some 
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Si sufficiently close to t. Hence m := inf { H{s, y) | s G [si,t), y G [x, 6+ (s)] } > 
by the continuity of H. Moreover, since V is continuous and V(t, y) = for all 
y £ )], it follows that 

\^iV{s,y)\<-[b+it-)-x] (3.32) 

for all s E [s2,t] and y S [x,6+(s)] where S2 S [si,t) is some value sufficiently close 
to t. Therefore for any s £ [s2,t) we find using (I3.22P and (|3.24p that 

V{s,x)= / / Vxx{s,z)dzdy (3.33) 

J X J y 

rb+{s) i-b+{s) 

= -2/ / {Vt + i^Vx + H){s,z)dzdy 

J X J y 

< 2^ ^ [^iV{s,y)-m{h+{s)-y)) dy 

< I {b+{t-)-x){b+is)-x) -m{b+is)-xf 

where in the second last inequality we used Vt > and in the last we used (j3.32p . 
Passing to the limit as s | t we see that V{t,x) < -f{b+{t-)-x) < which 
contradicts the fact that {t,x) £ D. We conclude that t ^+(i) is continuous on 
[0,T]. Note that this proof also shows that bj^{T) = = b^{T) since hj^{T) = = 
/i_(r) and V{x,T) = for all x £ ]R. 

8. Integral equations. We may now derive the integral equations ()3.5p and 
(13.6p + ()3.7p . From equation ()3.30p we see that 

\Vx{t,x)\ < 4{2Vt+\h\T) + 2E[Z] =:K<oo (3.34) 

for all {t,x) £ [0,T] xiR where Z := sup^gjj ^^,(-6^). Inserting this into ([3:22]) and 
using dSSD we obtain V^^ = -2{Vt+fiVx+H) < 2{\fi\ K-H) in C. If we let 

{l + \fi\K-H{t,z)) dzdy (3.35) 

for all {t, x) £ [0, T]xIR then we have Vxx < fxx on CU D° since \H\ < 1. Defining 
therefore the function F : [0, T] x IR ^ M hy F{t,x) = V{t,x) — f{t,x), we see 
that: (i) the map x F{t,x) is concave on each of the intervals (— oo,6_(t)), 
(6_(t),6+(t)) and (6+(t),oo) for every t £ [0,T]; (ii) the function F is C^'^ on 
C U -D°; (iii) the function Ft + fiF^ + ^Fxx is locally bounded on C U D°; and 
(iv) the maps t Fx{t,b-{t) ± ) = -fx{t,b-{t)) and t ^ Fx{t,b+{t) ± ) = 
—fx{t, b+{t)) are continuous on [0, T]. Since the boundaries 6_ and are monotone 
and consequently of bounded variation, we may apply the local time-space formula 
0] to F{t+s,B^_^J and Ito's formula to f{t+s,Bl^_^J since / is C^''^. Adding these 



t+r 
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two expressions and using (j3.22|) and (j3.24p we obtain 

Vit+s,B^^^) = V{t,x) (3.36) 

/(i?,%^{6_(t+r),6+(t+r)}) dr 

+ rv^{t + r,B^^^)l{Bj^^^.^{b.it+r),b+{t+r)}) dB, 
Jo 

+ lJ^{v,{t+r,B^^,+) - V,{t+r,B^^^ 

+ I J^{v,{t + r,B^_^,+) - K(t + r,i?,V-)^ 

= b+it + r)) dl\Xr{Bn 

= V{t,x) - [ H{t+r,B^^^.)l{b^{t+r) < B^^^ < 6+ (t+r)) dr + 
Jo 

under Pt,x for any (t,x) G [0,r]xiR and s G [0,T-t], where t-{B^') and i^+iBf") 
are the local times of B'^ on the curves 6_ and 6+ respectively, and Ms = Jq Vx{t+ 
r,B^j^^) dBt+r is a continuous martingale under Pt.x for s G [0,T — t]. Setting s = 
T — t, taking Et^x on both sides and using the optional sampling theorem we obtain 



/(i?,% = 6_(t+r)) di'f+ABn 



Vit,x) 



-t,x 



fT-t 

/ H{t + r,B>l^.,)l{b^{t+r)<B'^^,<b+{t + r)) dr 
Jo 



(3.37) 



which is exactly (13. 5h after interchanging the order of integration. Setting x equal 
to b-{t) and 6+(t) in ()3.37p we get 



T-t 



T-t 



-t,b-{t) 



-t,b+{t) 



H{t+r,B^^^.)l{b_{t+r)<B^^^.<b+{t + r)) dr = (3.38) 
H{t+r,B^_^^^)I{b-{t+r)<B^_^_^<b+{t+r))] dr = (3.39) 



which is exactly 



+ 



as claimed. 



9. Uniqueness. Lastly we show that 5_ and 6+ are the unique solutions to 
(|3.6p + (|3.7p in the class of continuous functions t b-{t) and t i— > 6+(t) on [0,T] 
satisfying 6_(t) < h_{t) and b+{t) > h+{t) for ah t in [0,r]. 

Take any two continuous functions c+and c_ on [0, T] which solve (j3.6p + (|3.7p and 
satisfy C-{t) < h-{t) and h+{t) < c+(t) for all t G [0, T]. Motivated by equation 
(13.371) above, define the continuous function : [0, T] x iR — > iR as 



U''{t,x) 



-t,x 



l-T-t 

/ H{t + r,B^!)l{c-{t + r) < B^^^. < c+(t + r)) dr 
Jo 



(3.40) 



and observe that c_ and c+ solving (j3.6p + (j3.7p means exactly that 
U^{t,c-{t)) = U^{t,c+it)) = for all t G [0,T]. Let := 
{ {t,x) G [0, T]xiR I X < c^{t) or X > c+(t)} so that Dc is closed and plays the 
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role of a 'stopping region' for c_ and C-^. To avoid confusion we will denote by Df, 
the original stopping set from (j3.3p defined by the functions 6_ and 
(i) We show that = on Dc- Since is Markov, the process 

A^, ■.= U%t+s,B^_^J+ r H{t+r,B^_^^) l{c-{t+r)<B^_^^<c+{t+r)) dr (3.41) 

^0 

is a martingale under Pf ^ for all s £ [0, T—t] and x £ M. Take any point (t, x) £ D^. 
and consider the stopping time 



c7c = inf{sG [0,T-t] I ^Z),} 



(3.42) 



under the measure P^^^,'- Since C/^ is zero on the curves c_|_ and c_ and [/^(T, x) = 
for all X £ M, we must have U'^{t + ac, B^_^^^J = 0. Inserting o"c in (j3.4ip . taking 
Pt^x expectations and using the optional sampling theorem (since H is bounded) 
we find that 



U'it,x) = Et,^ U'{t+ax,B^^^^ 







(3.43) 



showing that = on Dc as claimed. 

(ii) We show that C/'=(t,x) > V(t,x) for all {t,x) £ [0,T]x]R. To see this take 
any {t,x) £ [0,T] xJR and consider the stopping time 



rc = inf{s G [0,r-t] | B^^^ £ D^} 



(3.44) 



under Pt,x- We claim that U'^{t+Tc, B^_^^ ) = 0. Indeed, if (t, x) £ then Tc = so 
that U^{t,x) = by the argument above. Conversely if {t,x) ^ then the result 
follows since is zero on c-_ and c+ and U^{T,x) = for all x £ JR. Inserting Tc 
in (|3.4ip and using the optional sampling theorem, we see that 



U^{t,x) 



F(t+.,<J/(i3,%^Z),) ds 
F(t+s,i?f+J J >V{t,x) 



(3.45) 



where the second identity follows from the definition of Tc- We conclude that U'^ > 
y on [0, T] X iR as claimed. 

(iii) We show that Di^ C D^. Suppose this is not the case so that there exists 
some time t £ [0, T) for which either hj^{t) < c+(t) or b-{t) > c_(t). Choose any 
X > C-f-{t) OT X < c^{t) and consider the stopping time 



ab = inf{s£[0,T-t] | ^ A } 



(3.46) 



under the measure Pt,x- Replacing s with in (|3.36p and (|3.4ip and using the 
optional sampling theorem we find 



V{t + a,,B^^J =V{t,: 



U%t + ab,B'^ 



U%t,x) 



(3.47) 
(3.48) 



-t,x 



rcTb 

/ Hit+s,B^^,)l{B^^^^Dc) ds . 
Jo 
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Since {t, x) belongs to both and Dc it follows that U'^{t, x) = V{t, x) = 0, and 



the fact that W'it+ab, B^,^ ) > Vit+ab,B^, 



t+ah > 



implies 



-t,x 



rHit+s,B^^Jl{B^^^^D,) ds 
Jo 



< 0, 



(3.49) 



The assumption that either b^{t) < C4_(t) or b^{t) > c_(t) together with the 
continuity of the functions c_ , c+ , 6_ and 6+ means that there exists a small enough 
t < u <T so that 6+(s) < c+(s) or b-{s) > c_(s) for all s S [t, u]. Consequently the 
Pt,x probability of i?^ spending a strictly positive amount of time (w.r.t. Lebesgue 
measure) in either of these regions is strictly positive. Combined with the fact that 
both Db and Dc are contained in {H>0}, this forces the expectation above to be 
strictly positive and provides a contradiction. 

(iv) We show that Dc = D^. Suppose that this is not the case so that c+(t) < 
5+(i) or c_(t) > b-{t) for some t £ [0,T]. Choose any point x £ (c-|_(t), 6+(t)) or 
X G (6_(t), c_(i)) and consider the stopping time 



TD = mf{sG [0,T-t] I B^_^^ E A } 



(3.50) 



under Pt,x- Inserting td in (I3.36P and (I3.4ip . taking Pf^^; expectations and using 
the optional sampling theorem we obtain 



H{t+s,B^.Jds 



V{t,x) 



U'{t + TD,B^^^J =U'{t,x) 



-t,x 



H{t+s,B^^Jl{B^^,^Dc) ds 



(3.51) 
(3.52) 



Since D^ is contained in Dc and If^ = on Dc we must have W^it+To, = 0, 

and using the fact that If^ > V we find that 



-t,x 



rHit+s,B^^Jl{B^^,GDc) ds 
Jo 



< 0. 



(3.53) 



However as before the continuity of the boundaries b-,b^, c_ and c+ and the fact 
that Dc C {H > 0} forces the expectation to be strictly positive and provides a 
contradiction. We therefore conclude that c+(t) = b^{t) and c_(t) = b-{t) for all 

te[o,T]. □ 

Remark 1. On closer inspection of the proof of Lemma 1 it is clear that a 
similar argument could be made for a broad class of Markov processes. Indeed, in 
order to reduce the optimal prediction problem to an equivalent optimal stopping 
problem, the only difficulty lies in computing the conditional distribution of g. For 
any continuous time-homogenous Markov process X a simple application of the 
Markov property upon noting that maxt<,s<T Xg = ( maxo<s<T~i ^s) ° where 
Of is the shift operator shows that 



Pxig<t\J^t 



Px, i min Xs > O) if > 

^ 0<s<T-t ' 

if = 

Px, ( max Xs < O) if < 



(3.54) 
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for t G (0,T], and the problem reduces to computing the right-hand side of (j3.54p . 
Note that g < t in (j2.7p and (|2.8p can be replaced hy g < t so that the latter 
excludes the case Xf = iov t £ [0, T] as used in (j3.54p . 

For discontinuous Markov processes the problem is more involved as the repre- 
sentation ()3.54p above is no longer valid. In this case one can formulate an anal- 
ogous optimal prediction problem as follows. Let X = {Xt)o<t<T be a standard 
Markov process defined on a filtered probability space {J^t)o<t<T, Px) where 

Xq = X under P^^ for x € JR. Assume that is regular for itself in the sense that 
Po(ro = 0) = 1 where tq = inf { t G (0, T] | = } , let i = {£t)o<t<T denote the 
local time of X at (see P, p. 216]), and set 

5 = inf{tG [0,r] |4 = £t}. (3.55) 

Recalling that t i— > £f is continuous we see that ig = £t, and if X is continuous 
then g coincides with the last zero of X before T since t ^ it increases only when 
= for t G [0, r]. Inserting g from ()3.55p into (|2.2p we obtain an analogous 
optimal prediction problem which can be reduced to a standard optimal stopping 
problem. Indeed, since i is an additive functional it follows by the Markov property 
of X that 

PA9 <t\Tt) = P-A^t = ir I J't) = Px{£t = it + ir-t o et I Tt) (3.56) 
= PA^T^t o9t = 0\J^t) = PxA^T-t = 0) = F{T-t, Xt) 

where we set F{t,x) = Px{£t = ^) for (t, x) G [0,T] xJR. Having (|3.56p we can then 
proceed as in Lemma 1 above and define the value function V as in (j2.12p where 
H{t, x) = 2F{T-t, x)-l for (t, x) G [0, T]xlR. Note that t ^ F(t, x) is decreasing 
so that t ^ V{t, x) is decreasing and consequently if (t, x) belongs to the stopping 
set D then all points (i+s,x) belong to D for s G [0, T — t] with x £ M given and 
fixed. Note also that if X is continuous then 



F{T-t,x) = PAiT-t = 0) = { 



P^i min Xs>0) if x > 

^ 0<s<T-t ' 



max 

0<s<T-t 



if x = (3.57) 
X, < O) if X < 



for t G (0, T] and x G iR, which is in agreement with p.54p above. Finally, it is also 
evident that the previous considerations are by no means restricted to dimension 
one. 



Remark 2. Note that there are several options if one wishes to compute numeri- 
cal values for the value function V and the optimal stopping boundaries 6_ and 
The first is through the equations (j3.5p and (j3.6p + (j3.7p . which will involve the solu- 
tion of a coupled system of highly non-linear integral equations (this is the method 
applied to obtain Figure 1 above). Another is through the free-boundary problem 
(|3.22p - (j3.24p using a finite difference approach. One could also use a discrete-time 
discrete-space Markov chain approximation to and solve the analogous discrete 
problem using backward induction. The various methods each have their own pe- 
culiarities and strengths, and before trying to extract numbers from the problem, 
it is worth considering them all in order to choose the most appropriate one. 
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